Thermodynamics of a Hyperthermostable Carboxylesterase from Anoxybacillus geothermalis D9.

Hyperthermostable enzymes are very attractive biocatalysts due to their exceptional stability at extreme temperatures. Recently, a hyperthermostable carboxylesterase EstD9 from Anoxybacillus geothermalis D9 was biochemically characterized. The enzyme displayed remarkable stability at high temperature. In this study, we attempted to probe the conformational adaptability of EstD9 under extreme conditions via in silico approaches. Circular dichroism revealed that EstD9 generated new ß-sheets at 80 °C that make up the core of the hydrolase fold. Interestingly, the profiles of molecular dynamics simulation showed the lowest scores of radius of gyration and solvent accessible surface area (SASA) at 80 °C. Three loops were responsible for protecting the catalytic site, residing at the interface between the two domains. To further investigate molecular adaptation in extreme conditions, the intramolecular interactions of native structure were investigated. 18 hydrogen bond networks, 7 salt bridges, and 9 hydrophobic clusters were revealed within EstD9, which is higher than the reported thermostable carboxylesterase Est30. Collectively, the analysis indicates that intramolecular interactions and structural dynamics play distinct roles in preserving the overall EstD9 structure at elevated temperatures. This work is relevant to fundamental and applied research involving protein engineering of industrial thermostable enzymes.

Functional expression, purification, biochemical and biophysical characterizations, and molecular dynamics simulation of a histidine acid phosphatase from Saccharomyces cerevisiae.

A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110 kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6 s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75 min and 55 s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, ß-Turns, and ß-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.

Exploring the Antarctic aminopeptidase P from Pseudomonas sp. strain AMS3 through structural analysis and molecular dynamics simulation.

Aminopeptidase P (APPro) is a crucial metalloaminopeptidase involved in amino acid cleavage from peptide N-termini, playing essential roles as versatile biocatalysts with applications ranging from pharmaceuticals to industrial processes. Despite acknowledging its potential for catalysis in lower temperatures, detailed molecular basis and biotechnological implications in cold environments are yet to be explored. Therefore, this research aims to investigate the molecular mechanisms underlying the cold-adapted characteristics of APPro from Pseudomonas sp. strain AMS3 (AMS3-APPro) through a detailed analysis of its structure and dynamics. In this study, structure analysis and molecular dynamics (MD) simulation of a predicted model of AMS3-APPro has been performed at different temperatures to assess structural flexibility and thermostability across a temperature range of 0-60 °C over 100 ns. The MD simulation results revealed that the structure were able to remain stable at low temperatures. Increased temperatures present a potential threat to the overall stability of AMS3-APPro by disrupting the intricate hydrogen bond networks crucial for maintaining structural integrity, thereby increasing the likelihood of protein unfolding. While the metal binding site at the catalytic core exhibits resilience at higher temperatures, highlighting its local structural integrity, the overall enzyme structure undergoes fluctuations and potential denaturation. This extensive structural instability surpasses the localized stability observed at the metal binding site. Consequently, these assessments offer in-depth understanding of the cold-adapted characteristics of AMS3-APPro, highlighting its capability to uphold its native conformation and stability in low-temperature environments. In summary, this research provides valuable insights into the cold-adapted features of AMS3-APPro, suggesting its efficient operation in low thermal conditions, particularly relevant for potential biotechnological applications in cold environments.Communicated by Ramaswamy H. Sarma.

Effects of alcohol concentration and temperature on the dynamics and stability of mutant Staphylococcal lipase.

The stability and activity of lipase in organic media are important parameters in determining how quickly biocatalysis proceeds. This study aimed to examine the effects of two commonly used alcohols in industrial applications, methanol (MtOH) and ethanol (EtOH) on the conformational stability and catalytic activity of G210C lipase, a laboratory-evolved mutant of Staphylococcus epidermidis AT2 lipase. Simulation studies were performed using an open-form predicted structure under 30, 40 and 50% of MtOH and EtOH at 25 °C and 45 °C. The overall enzyme structure becomes more flexible with increasing concentration of MtOH and exhibited the highest flexibility in 40% EtOH. In EtOH, the movement of the lid was found to be temperature-dependent with a noticeable shift in the lid position at 45 °C. Lid opening was evidenced at 50% of MtOH and EtOH which was supported by the increase in SASA of hydrophobic residues of the lid and catalytic triad. The active site remained mostly intact. An open-closed lid transition was observed when the structure was re-simulated in water. Experimental evaluation of the lipase stability showed that the half-life reduced when the enzyme was treated with 40% (v/v) and 50% (v/v) of EtOH and MtOH respectively. The finding implies that a high concentration of alcohol and elevated temperature can induce the lid opening of lipase which could be essential for the activation of the enzyme, provided that the catalytic performance in the active site is not compromised.Communicated by Ramaswamy H. Sarma.

Enhancement in T1 lipase purification recovery using the novel construct pGEX4T1/His-T1.

The Geobacillus zalihae strain T1 produces a thermostable T1 lipase that could be used for industrial purposes. Previously, the GST-T1 lipase was purified through two chromatographic steps: affinity and ion exchange (IEX) but the recovery yield was only 33%. To improve the recovery yield to over 80%, the GST tag from the pGEX system was replaced with a poly-histidine at the N-terminal of the T1 lipase sequence. The novel construct of pGEX/His-T1 lipase was developed by site-directed mutagenesis, where the XbaI restriction site was introduced upstream of the GST tag, allowing the removal of tag via double digestion using XbaI and EcoRI (existing cutting site in the pGEX system). Fragment of 6 × His-T1 lipase fusion was synthesized, cloned into the pGEX4T1 system, and expressed in Escherichia coli BL21 (DE3) pLysS, resulting in lipase-specific activity at 236 U/mg. The single purification step of His-T1 lipase was successfully achieved using nickel Sepharose 6FF with an optimized concentration of 5 mM imidazole for binding, yielding the recovery of 98%, 1,353 U/mg lipase activity, and a 5.7-fold increase in purification fold. His-T1 lipase was characterized and was found to be stable at pH 5-9, active at 70 °C, and optimal at pH 9.

Carboxylic acid reductases: Structure, catalytic requirements, and applications in biotechnology.

Biocatalysts have been gaining extra attention in recent decades due to their industrial-relevance properties, which may hasten the transition to a cleaner environment. Carboxylic acid reductases (CARs) are large, multi-domain proteins that can catalyze the reduction of carboxylic acids to corresponding aldehydes, with the presence of adenosine triphosphate (ATP) and nicotinamide adenine dinucleotide phosphate (NADPH). This biocatalytic reaction is of great interest due to the abundance of carboxylic acids in nature and the ability of CAR to convert carboxylic acids to a wide range of aldehydes essentially needed as end products such as vanillin or reaction intermediates for several compounds production such as alcohols, alkanes, and amines. This modular enzyme, found in bacteria and fungi, demands an activation via post-translational modification by the phosphopantetheinyl transferase (PPTase). Recent advances in the characterization and structural studies of CARs revealed valuable information about the dynamics, mechanisms, and unique features of the enzymes. In this comprehensive review, we summarize the previous findings on the phylogeny, structural and mechanistic insight of the domains, post-translational modification requirement, strategies for the cofactors regeneration, the extensively broad aldehyde-related industrial application properties of CARs, as well as their recent immobilization approaches.

X-ray crystallography of mutant GDSL esterase S12A of Photobacterium marinum J15.

GDSL esterase is designated as a member of Family II of lipolytic enzymes known to catalyse the synthesis and hydrolysis of ester bonds. The enzyme possesses a highly conserved motif Ser-Gly-Asn-His in the four conserved blocks I, II, III and V respectively. The enzyme characteristics, such as region-, chemo-, and enantioselectivity, help in resolving the racemic mixture of single-isomer chiral drugs. Recently, crystal structure of GDSL esterase from Photobacterium J15 has been reported (PDB ID: 5XTU) but not in complex with substrate. Therefore, GDSL in complex with substrate could provide insights into the binding mode of substrate towards inactive form of GDSL esterase (S12A) and identify the hot spot residues for the designing of a better binding pocket. Insight into molecular mechanisms is limited due to the lack of crystal structure of GDSL esterase-substrate complex. In this paper, the crystallization of mutant GDSL esterase (S12A) (PDB ID: 8HWO) and its complex with butyric acid (PDB ID: 8HWP) are reported. The optimized structure would be vital in determining hot spot residue for GDSL esterase. This preliminary study provides an understanding of the interactions between enzymes and hydrolysed p-nitro-phenyl butyrate. The information could guide in the rational design of GDSL esterase in overcoming the medical limitations associated with racemic mixture.

Immobilization of Hyperthermostable Carboxylesterase EstD9 from Anoxybacillus geothermalis D9 onto Polymer Material and Its Physicochemical Properties.

Carboxylesterase has much to offer in the context of environmentally friendly and sustainable alternatives. However, due to the unstable properties of the enzyme in its free state, its application is severely limited. The present study aimed to immobilize hyperthermostable carboxylesterase from Anoxybacillus geothermalis D9 with improved stability and reusability. In this study, Seplite LX120 was chosen as the matrix for immobilizing EstD9 by adsorption. Fourier-transform infrared (FT-IR) spectroscopy verified the binding of EstD9 to the support. According to SEM imaging, the support surface was densely covered with the enzyme, indicating successful enzyme immobilization. BET analysis of the adsorption isotherm revealed reduction of the total surface area and pore volume of the Seplite LX120 after immobilization. The immobilized EstD9 showed broad thermal stability (10-100 °C) and pH tolerance (pH 6-9), with optimal temperature and pH of 80 °C and pH 7, respectively. Additionally, the immobilized EstD9 demonstrated improved stability towards a variety of 25% (v/v) organic solvents, with acetonitrile exhibiting the highest relative activity (281.04%). The bound enzyme exhibited better storage stability than the free enzyme, with more than 70% of residual activity being maintained over 11 weeks. Through immobilization, EstD9 can be reused for up to seven cycles. This study demonstrates the improvement of the operational stability and properties of the immobilized enzyme for better practical applications.

In silico structural exploration of serine protease from a CTG-clade yeast Meyerozyma guilliermondii strain SO.

In eukaryotes, serine proteases are cellular localized hydrolases reported to regulate essential biological reactions. Improved industrial applications of proteins are aided by prediction and analysis of their 3-dimensional structures (3D). A serine protease was identified from CTG-clade yeast Meyerozyma guilliermondii strain SO and its 3D structure as well as its catalytic attributes have not been fully understood yet, thus we seek to report on the catalytic mechanism of M. guilliermondii strain SO MgPRB1 using substrate PMSF via in silico docking as well as its stability by way of disulfide bonds formation. Herein, bioinformatics tools and techniques were used to predict, validate and analyze the possible changes of CUG ambiguity (if any) in strain SO using template PDB ID: 3F7O. Structural assessments confirmed the classic catalytic triad Asp305, His337, and Ser499. Superimposition of MgPRB1 and template 3F7O structures revealed the unlinked cysteine residues between Cys341, Cys440, Cys471 and Cys506 of MgPRB1 compared to template 3F7O with two disulfide bonds formation, which confers structural stability. In conclusion, serine protease structure from strain SO was successfully predicted and studies towards understanding at the molecular level may be undertaken for its potential applications in the degradation of peptide bonds.

Knotting terminal ends of mutant T1 lipase with disulfide bond improved structure rigidity and stability.

Lipase biocatalysts offer unique properties which are often impaired by low thermal and methanol stability. In this study, the rational design was employed to engineer a disulfide bond in the protein structure of Geobacillus zalihae T1 lipase in order to improve its stability. The selection of targeted disulfide bond sites was based on analysis of protein spatial configuration and change of Gibbs free energy. Two mutation points (S2C and A384C) were generated to rigidify the N-terminal and C-terminal regions of T1 lipase. The results showed the mutant 2DC lipase improved methanol stability from 35 to 40% (v/v) after 30 min of pre-incubation. Enhancement in thermostability for the mutant 2DC lipase at 70 °C and 75 °C showed higher half-life at 70 °C and 75 °C for 30 min and 52 min, respectively. The mutant 2DC lipase maintained the same optimum temperature (70 °C) as T1 lipase, while thermally induced unfolding showed the mutant maintained higher rigidity. The kcat/Km values demonstrated a relatively small difference between the T1 lipase (WT) and 2DC lipase (mutant). The kcat/Km (s-1 mM-1) of the T1 and 2DC showed values of 13,043 ± 224 and 13,047 ± 312, respectively. X-ray diffraction of 2DC lipase crystal structure with a resolution of 2.04 Å revealed that the introduced single disulfide bond did not lower initial structural interactions within the residues. Enhanced methanol and thermal stability are suggested to be strongly related to the newly disulfide bridge formation and the enhanced compactness and rigidity of the mutant structure. KEY POINTS: • Protein engineering via rational design revealed relative improved enzymatic performance. • The presence of disulfide bond impacts on the rigidity and structural function of proteins. • X-ray crystallography reveals structural changes accompanying protein modification.

Recent advances in simultaneous thermostability-activity improvement of industrial enzymes through structure modification.

Engineered thermostable microbial enzymes are widely employed to catalyze chemical reactions in numerous industrial sectors. Although high thermostability is a prerequisite of industrial applications, enzyme activity is usually sacrificed during thermostability improvement. Therefore, it is vital to select the common and compatible strategies between thermostability and activity improvement to reduce mutants̕ libraries and screening time. Three functional protein engineering approaches, including directed evolution, rational design, and semi-rational design, are employed to manipulate protein structure on a genetic basis. From a structural standpoint, integrative strategies such as increasing substrate affinity; introducing electrostatic interaction; removing steric hindrance; increasing flexibility of the active site; N- and C-terminal engineering; and increasing intramolecular and intermolecular hydrophobic interactions are well-known to improve simultaneous activity and thermostability. The current review aims to analyze relevant strategies to improve thermostability and activity simultaneously to circumvent the thermostability and activity trade-off of industrial enzymes.

Bacteria consortia enhanced hydrocarbon degradation of waxy crude oil.

Waxy crude oil is a problem to the oil and gas industry because wax deposition in pipelines reduces the quality of the crude oil. Currently, the industry uses chemicals to solve the problem but it is not environmentally friendly. As an alternative, the biodegradation approach is one of the options. Previously eleven thermophilic bacteria were isolated and exhibited high ability to degrade hydrocarbon up to 70% of waxy crude oil. However, despite the successful study on these single bacteria strains, it is believed that biodegradation of paraffin wax requires more than a single species. Five consortia were developed based on the biodegradation efficiency of 11 bacterial strains. Consortium 3 showed the highest biodegradation (77.77%) with more long-chain alkane degraded throughout the incubation compared to other consortia. Enhancement of hydrocarbon degradation was observed for all consortia especially in long chain alkane (C18-C40). Consortium 3 exhibited higher alkane monooxygenase, alcohol dehydrogenase, lipase, and esterase activities. Moreover, the dominant bacteria in the consortia were determined by denaturing gradient gel electrophoresis (DGGE), which showed the domination of genera Geobacillus, Parageobacillus, and Anoxybacillus. It can be concluded that the bacterial consortia showed higher biodegradation and improved degrading more long-chain hydrocarbon compared to a single isolate.

Characterization of Carboxylic Acid Reductase from Mycobacterium phlei Immobilized onto Seplite LX120.

A multi-domain oxidoreductase, carboxylic acid reductase (CAR), can catalyze the one-step reduction of carboxylic acid to aldehyde. This study aimed to immobilize bacterial CAR from a moderate thermophile Mycobacterium phlei (MpCAR). It was the first work reported on immobilizing bacterial CAR onto a polymeric support, Seplite LX120, via simple adsorption. Immobilization time and protein load were optimized for MpCAR immobilization. The immobilized MpCAR showed optimal activity at 60 °C and pH 9. It was stable over a wide range of temperatures (10 to 100 °C) and pHs (4-11), retaining more than 50% of its activity. The immobilized MpCAR also showed stability in polar solvents. The adsorption of MpCAR onto the support was confirmed by Scanning Electron Microscopy (SEM), Fourier-Transform Infrared (FTIR) spectroscopy, and Brunauer-Emmett-Teller (BET) analysis. The immobilized MpCAR could be stored for up to 6 weeks at 4 °C and 3 weeks at 25 °C. Immobilized MpCAR showed great operational stability, as 59.68% of its activity was preserved after 10 assay cycles. The immobilized MpCAR could also convert approximately 2.6 mM of benzoic acid to benzaldehyde at 60 °C. The successfully immobilized MpCAR on Seplite LX120 exhibited improved properties that benefit green industrial processes.

A new hyper-thermostable carboxylesterase from Anoxybacillus geothermalis D9.

Carboxylesterases are attractive biocatalysts for various industrial applications, especially hyperthermophilic carboxylesterases, due to their high tolerance toward extreme environments. Such ability confers many advantages, including cost-effectiveness and an increased manufacturing rate. In the current work, we first described the characterization of EstD9, a new carboxylesterase from thermophilic Anoxybacillus geothermalis D9. Sequence analysis of EstD9 revealed a significant identity (80 %) with thermophilic Est30 and a catalytic triad, composed of Ser93-His22-Asp193. As the protein sequence contained a conserved pentapeptide (GLSLG), EstD9 could be proposed as a new member of family XIII. The putative carboxylesterase was recombinantly expressed in E. coli BL21 (DE3) with a molecular mass of 28 kDa and successfully purified via affinity chromatography with recovery of 88.36 %. Using p-nitrophenyl butyrate, EstD9 presented excellent stability at high temperature range (70 °C-100 °C) and a broad pH tolerance (pH 6-9), with optimal activity at 80 °C and pH 7. Notably, EstD9 activity was stimulated in the presence of 1-propanol and DMSO with 107.8 % and 108.9 % relative activities, respectively. The purified EstD9 maintained 60 % residual activity after 30 min exposure to various surfactants and metal ions. Additionally, the inhibition studies demonstrated strong deactivation by phenylmethylsulfonyl fluoride, dithiothreitol, and ß-mercaptoethanol. The estimated Tm value was 72.12 °C. Unlike typical carboxylesterases, in silico 3D model of EstD9 disclosed a topological α/ß hydrolase fold with a small α-helix cap. The enzymatic properties of EstD9 suggest this enzyme to be a highly suitable catalyst for industrial bioprocesses under harsh conditions.

Membrane fatty acid desaturase: biosynthesis, mechanism, and architecture.

Fatty acid desaturase catalyzes the desaturation reactions by inserting double bonds into the fatty acyl chain, producing unsaturated fatty acids, which play a vital part in the synthesis of polyunsaturated fatty acids. Though soluble fatty acid desaturases have been described extensively in advanced organisms, there are very limited studies of membrane fatty acid desaturases due to their difficulties in producing a sufficient amount of recombinant desaturases. However, the advancement of technology has shown substantial progress towards the development of elucidating crystal structures of membrane fatty acid desaturase, thus, allowing modification of structure to be manipulated. Understanding the structure, mechanism, and biosynthesis of fatty acid desaturase lay a foundation for the potential production of various strategies associated with alteration and modifications of polyunsaturated fatty acids. This manuscript presents the current state of knowledge and understanding about the structure, mechanisms, and biosynthesis of fatty acid desaturase. In addition, the role of unsaturated fatty acid desaturases in health and diseases is also encompassed. This will be useful in understanding the molecular basis and structural protein of fatty acid desaturase that are significant for the advancement of therapeutic strategies associated with the improvement of health status. KEY POINTS: • Current state of knowledge and understanding about the biosynthesis, mechanisms, and structure of fatty acid desaturase. • The role of unsaturated fatty acid desaturase. • The molecular basis and structural protein elucidated the crystal structure of fatty acid desaturase.

Structural interpretations of a flexible cold-active AMS8 lipase by combining small-angle X-ray scattering and molecular dynamics simulation (SAXS-MD).

Determining structure of highly flexible protein with multiple conformations can be challenging. This paper aims to combine molecular dynamics (MD) and small angle X-ray diffraction (SAX) techniques as a solution to overcome issues related to protein conformation in hardly crystallized protein. Based on prior studies, a cold-active lipase AMS8 was simulated in solvents showing stability in its N-terminal and high flexibility in its C-terminal. However, MD in its own algorithm could not explain the basis of macromolecule conformational transitions or changes related to protein through folding. Hence, by combining SAXS with MD, it is possible to understand the structure of flexible AMS8 lipase in natural space. Based on the findings, SAXS ab-initio model of AMS8 lipase was identified as a monomeric protein in which the optimized model of cold-active lipase AMS8 derived from SAXS data was found to be aligned with AMS8 homology model under series of MD timeframe.

Structure Prediction and Characterization of Thermostable Aldehyde Dehydrogenase from Newly Isolated Anoxybacillus geothermalis Strain D9.

In nature, aldehyde dehydrogenase (ALDH) is widely distributed and mainly involved in the oxidation of aldehydes. Thermostability is one of the key features for industrial enzymes. The ability of enzymes to withstand a high operating temperature offers many advantages, including enhancing productivity in industries. This study was conducted to understand the structural and biochemical features of ALDH from thermophilic bacterium, Anoxybacillus geothermalis strain D9. The 3D structure of A. geothermalis ALDH was predicted by YASARA software and composed of 24.3% ß-sheet located at the center core region. The gene, which encodes 504 amino acids with a molecular weight of ~56 kDa, was cloned into pET51b(+) and expressed in E.coli Transetta (DE3). The purified A. geothermalis ALDH showed remarkable thermostability with optimum temperature at 60 °C and stable at 70 °C for 1 h. The melting point of the A. geothermalis ALDH is at 65.9 °C. Metal ions such as Fe3+ ions inhibited the enzyme activity, while Li+ and Mg2+ enhanced by 38.83% and 105.83%, respectively. Additionally, this enzyme showed tolerance to most non-polar organic solvents tested (xylene, n-dedocane, n-tetradecane, n-hexadecane) in a concentration of 25% v/v. These findings have generally improved the understanding of thermostable A. geothermalis ALDH so it can be widely used in the industry.

Thermostability engineering of industrial enzymes through structure modification.

Thermostability is an essential requirement of enzymes in the industrial processes to catalyze the reactions at high temperatures; thus, enzyme engineering through directed evolution, semi-rational design and rational design are commonly employed to construct desired thermostable mutants. Several strategies are implemented to fulfill enzymes' thermostability demand including decreasing the entropy of the unfolded state through substitutions Gly → Xxx or Xxx → Pro, hydrogen bond, salt bridge, introducing two different simultaneous interactions through single mutant, hydrophobic interaction, filling the hydrophobic cavity core, decreasing surface hydrophobicity, truncating loop, aromatic-aromatic interaction and introducing positively charged residues to enzyme surface. In the current review, horizons about compatibility between secondary structures and substitutions at preferable structural positions to generate the most desirable thermostability in industrial enzymes are broadened. KEY POINTS: • Protein engineering is a powerful tool for generating thermostable industrial enzymes. • Directed evolution and rational design are practical approaches in enzyme engineering. • Substitutions in preferable structural positions can increase thermostability.

Versatility of subtilisin: A review on structure, characteristics, and applications.

Due to its thermostability and high pH compatibility, subtilisin is most known for its role as an additive for detergents in which it is categorized as a serine protease according to MEROPS database. Subtilisin is typically isolated from various bacterial species of the Bacillus genus such as Bacillus subtilis, B. amyloliquefaciens, B. licheniformis, and various other organisms. It is composed of 268-275 amino acid residues and is initially secreted in the precursor form, preprosubtilisin, which is composed of 29-residues signal peptide, 77-residues propeptide, and 275-residues active subtilisin. Subtilisin is known for the presence of high and low affinity calcium binding sites in its structure. Native subtilisin has general properties of thermostability, tolerance to neutral to high pH, broad specificity, and calcium-dependent stability, which contribute to the versatility of subtilisin applicability. Through protein engineering and immobilization technologies, many variants of subtilisin have been generated, which increase the applicability of subtilisin in various industries including detergent, food processing and packaging, synthesis of inhibitory peptides, therapeutic, and waste management applications.